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ABSTRACT 

The structure of the broad emission line region (BLR) in active galactic nuclei (AGN) 
remains unclear. We test in this paper a flattened configuration model for BLR. The 
virial theorem, by taking into account the disc shape of BLR, allows us to get a 
direct connection between the mass of a supermassive black hole (SMBH) and the 
inclination angle of the accretion flow. The inclination angle itself is derived from the 
spectropolarimetric data on broad emission lines using the theory for the generation of 
polarized radiation developed by Sobolev and Chandrasekhar. As the result, the new 
estimates of SMBH masses in AGN with measured polarization of BLR are presented. 
It is crucial that the polarimetric data allow also to determine the value of the virial 
coefficient that is essential for determining SMBH masses. 

Key words: polarization, accretion disc, supermassive black holes, active galactic 
nuclei. 


1 INTRODUCTION 

It is generally accepted that AGNs are powered by accretion 
to a SMBH and the broad emission lines seen in type I AGN 
are produced as the result of photoionization of gas in the 
BLR. However the structure and kinematics of the BLR is 
a key problem in modern astrophysics. Broad emission lines 
are emitted in the vicinity of SMBH in AGN, but this re¬ 
gion is not resolved in observations. Properties of the broad 
emission lines are used to estimate the mass of the central 
SMBH. The traditional method for estimating the SMBH 
mass is to use the virial theorem. It means that the mass 
of the SMBH can be e stimated as ijVestergaard fc Peterson! 
l2006l:lFine et al.ll2008l) : 

= ( 1 ) 

where Mbh is the mass of a black hole, / is a virial param¬ 
eter that defines the geometry, velocity field and orientation 
of BLR, Rblr is the radius of BLR and Vblr is the velocity 
dispersion that is measured as the full width of the emis¬ 
sion line in the variable spectrum. The BLR radius Rblr 
is usually given by the time delay between continuum and 
emission line variations. 

There are var i ous ap proaches for determining the value 
of /. iLabita et al.l l|200fil l found that in quasars an isotropic 
BLR fails to reproduce the observed line widths and shapes, 
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and a disk model is preferred. A disk like geo metry for the 
BLR has been proposed by several authors ijPecarli et al.l 
I 2 OO 8 I) . Some autho rs suggested that t he BLR can not be 
completely flat (e.g. lCollin et ah] (l2006l )). It means that disk 
may have a finite half thickness H, or a profile with H in¬ 
creasing more than linearly with the disk radius. Other mod¬ 
els pr opose the existence of warped disks (iTremaine fc Da\^ 

I 2 OI 4 I) . 

According to iKo llatschnv. Zetzl fc Dietrichl ll2006l) : 
iKollatschnv fc Zetzll 1 20131 ) the hydrogen lines are emitted 
in a more flattened configuration in comparison to the highly 
ionized lines. Hp lines originate at heights of 0.7 to 1.6 
light days and at distances of 1.4 to 2.4 light days with 
height /distance ratios H/R ~ 0.07 — 0.5. iPancoast et al.l 
ll2013t) found that the geometry of the BLR is generally 
a thic k disk viewed close to face on. lEracleous fc HalpernI 
l|l994h have f ound that the inclination a ngle of BLR is 
24° —30°, and lEracleous. Halpern fc Livid ll 19961) suggested 
that the inclination of the BLR can be i = 19° — 42°. 

The virial coefficient / depends strongly on the BLR 
geometry, veloci ty field and orientati o n. Ma ny authors used 


the value f ~ 1. Peterson &: Wand^ ll 19991 ) found / = 3/4. 


lOnken et al.l ll2004l) found the mean value of the virial coef¬ 
ficient / = 14. iMcLure fc Dunlonl 1 I 2 OOII) have shown that 
for a disk inclined at an angle i to the observer the virial 
coefficient value is 

In this paper we adopt the disk like model for the BLR 
of Seyfert galaxies and the expression for / given by ([2]). The 
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value of the inclination angle can be determined from po- 
larimetric obser vations using the standard Chandr asekhar- 
Sobolev theory (|Chandrasekhaill 195(11 : ISobolevIllQh.lI l of mul¬ 
tiple scatterings of the radiation on free electrons and 
Rayleigh scattering on gas molecules and small dust par¬ 
ticles. According to these classical works, the polarization 
degree of scattered radiation depends strongly on the in¬ 
clination angle. The scattered radiation has the maximum 
linear polarization Pi = 11.7% when the line of sight is per¬ 
pendicular to th e normal to the semi -infin ite atmosphere 
(Milne probleml. IChandrasekharl lll950h and lSobolevI (Il963l l 
presented the solution of the so-called Milne problem, i.e. 
multiple scattering of light in optically thick flattened at¬ 
mospheres. The Milne problem corresponds to the propaga¬ 
tion and scattering of light in optically thick disk-like region, 
i.e. this solution can be applied to BLR. The idea of infer¬ 
ring the inclination of a black hole accretion disk from ob- 
servation of its polarized continu um has been suggested by 
iLi, Naravan fc McClintod3 (l2009l l and developed for deter¬ 
mining the orientation of the X-ray producing inner region 
of the accretion disk around a black hole in X-ray binary sys¬ 
tems. Another id ea for using spectropol arimetric data has 
been suggested by I Afanasiev et ahl ([201^. The authors used 
the spectropolarimetric observation of the line profile for Sy 
1.5 galaxy Mrk 6 and assumed that the scattering region 
is loca ted in the inner part of the torus. Recently iMarinl 
Il20l4l . based on archival data, reported the first compila¬ 
tion of 53 estimated AGN inclinations matched with ultra¬ 
violet/optical continuum polarization measurements. 

But these data are obtained in the widebands of wave¬ 
length and do not include the data on the broad Ha and Hp 
emission that are used for determination of the masses of 
SMBH. The values of the inclination angles are obtained in 
many cases by indirect method s. For example, IWu fc HanI 
ll200lll and IZhang fc Wul (l2002l l estimated the orientation 
angles i assuming a definite mass/velocity dispersion be¬ 
tween AGN and regular galaxies where AGN are located. 
For some Seyfert galaxies the inclination estimations are ob¬ 
tained from optical polarization data that are produced in 
contin uum and can not b e direc tly related to BLR. Accord¬ 
ing to lTremaine fc DavisI (l2014l 'l the warped disks can exist 
in AGN and it means that BLR and the standard accretion 
disk, located closely to SMBH, can have different inclina¬ 
tions. This phenomenon is the reason why we used the po- 
larime tric data of Ha emission line presented bv ISmith et al.l 
(1200211 . 

We u se the theory of mult i ple scattering s of polarized 
radiation (iGhandrasekharll 1950l : ISobolevlll963l l and the disk 
like model for the BLR. As the result we estimate the values 
of the virial factor and the mass values for SMBH in AGN. 
W e have in our order the detailed atlas of the values of the 
polarization degree and its dependence on the inclination 
angle i with the step of 0.005 for /r = cosi. Besides it is 
convenient in some cases to use the an alytical formula for th e 
polarization degree Pi{n) obtained in ISilant’ev et all ll2010ll : 


oz 1 - 2.2399At -f 2.9377^^ - 1.6978 ai® 

1 + 2.2001^ + 0.062p2 - 0.1988/r3 ' 


(3) 


The determination of the SMBH mass depends strongly 
on the virial coefficient. According to Eq.(I21) the virial co- 
ef ficient depends on the incli nation angle. T he calculations 
of iGhandrasekharl lll950h and ISobolevI J 19631 ') of the degree 


of polarization of multiple scattered radiation coming out 
from the flattened atmosphere is strongly dependent on the 
inclination angle i. Thus, according to Eq. © the determi¬ 
nation from polarization data the value of the inclination 
angle allows us to determine the virial coefficient /. 


2 DETERMINING THE VIRIAL 

COEFFICIENT 

For determining the virial coefficient we use Eq. ([2J . The po- 
larimetric data that are necessary for determining the virial 
coefficient a r e pre sented in the spectropolarimetric atlas of 
ISmith et al.1 (l2002f) . They obtained the values of polarization 
degree and position angles for 36 Type I Seyfert galaxies dur¬ 
ing a number of different runs at the Anglo-Australian and 
William H erschel telescopes. From 36 objects presented in 
the atlas of ISmith et al.1 (l2002l l approximately 13 AGNs have 
the equal polarization degree for Ha line and continuum and 
approximately equal values of the position angle within the 
error limits. For most of the observed objects there is dif¬ 
ference between the values of the polarization degree and 
the position angle. This fact can testify to the difference in 
the inclinations between disk shaped BLR and the accre- 
tion disk that can be desc ribed by Shakura-Sunyaev model 
llShakura fc Sunvaevlll97^ L We shall show below that this 
conclusion is confirmed by polarimetric observations. Con¬ 
tinuum polarization degree of AGN from Palomar-Green 
catalogue have been measu red on 6-m telescope of Special 
Astrophysical Observatory (lAfanasiev et al.ll201ll i. 

Our results for the virial coefficient / are presented at 
Table 1. In the second column of this table we present the 
values f or the BLR i n clinat ions obtained from polarimetric 
data of ISmith et al.l (l2002l l . The virial factor for most of 
the objects from t he Table 1 is / 1 with in error limits, 

in accordance with iMcLure fc Jarvi^ (l2002|j . For some ob¬ 
jects, for example, Mrk 841, Mrk 896, Mrk 926, NGC 3516, 
NGC 3783, NGC 5548 the virial coefficient is in better 
agree ment with the mean value of / = 1.4 dOnken et al.l 
I 2 OO 4 II . For a number of objects including Fairall 51, Mrk 6, 
MC 1849.2-78.32, NGC 6814, UGC 3478 and WAS 45 the 
inclination angle exceeds considerably the value i = 30°, and 
these objects have values of the virial coefficient of / < 1.0. 
This fact can mean that fo r BLR the model of random orbits 
dPeterson fc Wanddlll999ll is more realistic. 


3 DETERMINING SMBH MASSES 


After deriving the virial coefficient, it is possible to deter¬ 
mine the SMBH mass in these particular AGN using Eq.((T]). 
For a disk shape BLR, Eq. m takes a form: 


where F WHM is the observed full width of the emis¬ 
sion line dVestergaard fc PetersonI 20061: Wang et alJ l2009l : 
IHo. Darling fc Greened 20081 : Feng. Shen fc Lil 2014ll . and 
Rg = GMbh/P is the gravitational radius. The quan¬ 
tity RblrVblr/G is called the ’’virial product”(VP). This 
quantity is based on the two observable quantities: BLR ra¬ 
dius and emission fine width and has units of mass. The 
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Table 1. The values of the virial coefficient / determined via 
measured pola rization of B r oad H ex, emission of AGN. Polariza¬ 
tion data from ISmith et aT] 1I2OO2I) . i is the inclination angle. 


Object 

i [deg] 

/ 

Akn 120 

29 ±1 

i.UD4_Q Q34 

Akn 564 

26tl 

9c:+0.16 
-L-^O_0 17 

ESO 012-021 

34 ±4 


Fairall 51 

77.5 ±0.5 

0.262 ±0.001 

I Zw 1 

26 ± 1 

1 Qi r + 0.095 

T. o±o_Q 085 

ESO 141-G35 

43.5t?;o 

r» e:o+0.032 

0»-00_o 015 

KUV 18217-1-6419 

16+1 

0 ig+0.45 

Mrk 6 

41+0.5 

'J-Oo_o 01 

Mrk 279 

2lti 

2 0*^^ ® 

^•'^-0.4 

Mrk 290 

29 ±2 

-*-•^'-0.11 

Mrk 304 

37-5lns 

0.6651°;°®! 

Mrk 335 

09 O+0-7 
3^.o_q g 

0.86l°;°t 

Mrk 509 

31.0++° 

0 94'*"® *^^ 
^•^^_0.04 

Mrk 705 

25.0l°;« 

1 r)K+0.33 
l.UO_o 19 

Mrk 841 


-| qq+0.29 
1.3(5_o 21 

Mrk 871 

36.512-5 

n 7+0.16 

^•'-0.07 

Mrk 876 

28.01+® 

1 1 7+0.21 
-*-•-*^'-0.14 

Mrk 896 

2215 

1 ^4+1.38 
'^-0.51 

Mrk 915 

30.oil]® 

1 rv + 0.02 

-^•^-0.19 

Mrk 926 

21.0l°;° 

1 95 

i.yo.u_o 40 

Mrk 985 

35.0l+° 

n 75+0-04 
'J-' O_o 02 

MS 1849.2-78.32 

50.5tl:l 

0 42+0-02 
'J-^^_0.02 

NGC 3516 

18 ±2.5 

^•0-0.45 

NGC 3783 

22.0ll° 

1 79+0-26 
'^-0.1 

NGC 4051 

31 ±3 

0 96+^'^'^ 
u.yo_o 15 

NGC 4593 

33.012;® 

0.80 ±0.06 

NGC 5548 

21.o15;® 

i.y(5_0 13 

NGC 6104 

30.0l°;® 

i.uu_o 18 

NGC 6814 

54 ± 1 

^ qq+u.u.u5 
U.3o_o 10 

NGC 7213 

21"*"'^ 

1.95 ±0.8 

NGC 7469 

isli 

0 70+u.Dy 

-J- ''5-0.82 

NGC 7603 

23.o12;° 

1 63+0-28 
±.uo_o 22 

PG 1211-1-143 

15 ±4 

3 62+’^' ^ 
o.oz_i 28 

UGC 3478 

45.5 ± 1.5 

0.49 ±0.3 

WAS 45 

61.512;° 

0.32 ±0.01 


dimensionless factor / in Eq. o depends on the structure, 
velocity field and inclination of BLR and it is different for 
each AGN. For the flattened, disk-like structure of BLR size 
is measured usually from reverberation mapping via the time 
lag between the broad line emission and continuum variabil¬ 
ities. Vblr is measured usually as full width at half maxi¬ 
mum (FWHM) or as the line width that is usually charac¬ 
terized by the broad line dispersion. The full width of half 
maximum (FWHM) is used in many papers. The value of 
sin i can be derived from the data of the degree of polariza¬ 
tion that is strongly dependent on the inclination angle. We 
have the results of the detailed numerical calculations of the 
degree of polarization Pi (p) for the radiation scattered in the 
optically thick plane parallel atmosphere. These calculations 
are made in the framework of the classical Chandrasekhar- 


Table 2. The masses of SMBH in AGNs determined via measured 
polarization of broad Ha emission. 


Object 

Mbh/Mq 

Mbh/Mq 


(observations) 

(literature) 

Akn 1202 

(3.23l°;!!) X 10® 

(4.49 ± 0.93) X 10® 

Akn 564! 

(I.2IH9) X 10® 

~ 1.1 X 10° 

Fairall 9^ 

(1.34l°+«) X 10® 

(2.55 ±0.56) X 10® 

I Zw l2 

Iq7.6±0.17 

7 441+°'°°^ 

10 -0.119 

Mrk 6° 

(1.09l° °!) X 10® 

(1.36 ±0.12) X 10® 

Mrk 279^ 

(8.13+24) X 10! 

(15.2+2|) X 10! 

Mrk 290° 

(3.94 ±0.19) X 10’! 

(2.43 ± 0.37) X 10'! 

Mrk 304^ 

10^'^-0.02 

8 511+0 093 

10».Oii_o 113 

Mrk 335! 

(1.56l+») X 10! 

(1.42 ±0.37) X 10'! 

Mrk 509° 

(1.35 ± 0.12) X 10® 

(1.39 ±0.12) X 10® 

Mrk 705^! 

7 07+°-“ 

10 -0.09 

]^q6.79+0.5 

Mrk 8412 

108.55+0.1 

Q e/oq + 0.079 

10 -0.052 

Mrk 871"! 

7 04+°-°° 

10 -0.06 

]^q7.08+0.5 

Mrk 8762 

^ ,-7+0.19 
10^-*^'-0.52 

q 1 oq+0.096 

10 -0.122 

Mrk 896® 

]^q7.07±0.06 

io+°i 

Mrk 926° 

Q O+0.19 

10^-S-o.ii 

108.36+0.02 

Mrk 985’! 

3.18 X Uf 

5.71 X 10’! 

NGC 3516® 

]^q8.06±0.27 

7 q^+0.04 

107.88_o 03 

NGC 3783° 

]^q7.7±0.11 

7 47+0-07 

10''^'-0.09 

NGC 4051!° 

(1.64t+!) X 10° 

(1-58+;®°) X 10° 

NGC 4593!° 

(8.25+!°) X 10° 

(9.8 ±2.1) X 10° 

NGC 55482 

(7.84t+°) X 10! 

(7.827 ±0.017) X Vf 

NGC 6IO4!! 

10+!°-“o:“o^8 

10+°9 

NGC 6814'* 

6 94+0-077 

10O.y4i_0 09 

]^q7.02+0.5 

NGC 7213'* 

f\ OO+0.84 

106.83_o 33 

106.88+0.5 

NGC 7469*° 

7 54+0.17 
10'-^^-0.22 

]^q7.19+0.13 

PG 1211±1432 

Q 04 + 0.29 

108.34_o 20 

7 Q61 +0.082 

10f.yDi_o 101 


(11 iRevnoldj ll2013l l: I2l I Vestergaard &: PetersonI (120061) : 

('3'l lFeng. Shen &: hi (120141 '): ('4 ') IHo. Darling fc Green j 1I2OO8II : 

(5) IShankar et al.l ll2012l i: (6) IWinter et al.l ll20inl '): 

171 IZiolkowskil 1I2OO8II : ISl IWu fc HanI ll200lh : 

(9)i rennema 3 ll201.'jl : (10) IWang et al.l ll200flll : 

(11) IWoo et ni ll20l4 l: 

Sobolev theory l|Chandrasekhail Il95fll : ISobolevI Il963h with 
the step for p being equal to 0.005. 

The results of our calculations of the SMBH masses are 
presented in Table 2. In the last column of the table the pub¬ 
lished values for the SMBH mass are presented. Our values 
show a good overall agreement with the previous determi¬ 
nations of black hole masses, obtained by various methods. 
For a number of objects there is a difference. For Fairall 9 
our estimate of the SMBH mass looks lower, but our up¬ 
per limit value co incides wi t h the low-order limit value of 
BH, presented bv iRevnoldsI ll2ni3h . For Ark 120 the situa¬ 
tion looks as the same and there is the coincidence of our 
estimated value with the low-order value of BH mass, pre- 
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sented by IVestergaard fc Peterson! (l2006l l. The similar sit¬ 
uation occurs also for NGC 3516, NGC 3783, NGC 5548 
and Mrk 279. Obly for Mrk 290 our low-order estimate 1.34 
times higher than the u pper limit of BH mass, presented by 
iFeng. Shen fc Lil ll2014h . 

The difference between our and previous estimates of 
mass for some black holes can be associated with the 
real determination of the virial parameter / from Eq. O- 
iMcLure fc Dunlop ll200lll value of / (Eq.(l2])) allows us to 
determine directly this parameter but only in the situation 
when the inclination of the BLR is certain. Polarimetric ob¬ 
servations have preference because the value of the polariza¬ 
tion degree is directly associated with the inclination angle 
value, especially, for the standard Chandrasekhar-Sobolev 
theory of the generation of polarization in the plane-parallel 
atmosphere. Unfortunately, other methods determining in¬ 
clination angle of the BLR and accretion disk are consider¬ 
ably uncertain. For example, it is used determination of the 
accretion disk inclination to the line of sight for a sample 
of AGNs from their bulge stellar velocity dispersion, based 
on suggesting a Keplerian mass/vel ocity di s persio n between 
AGN and host galaxies (details in iMari 3 (|2014lU . Namely 


the use of the bulge stellar velocity dispersion provides the 
difference in the real estimate of BH masses. In many cases 
the BH mass is estimated suggesting that the virial coeffi¬ 
cient / = 1. 

Another problem associated with the determination of 
the BLR size. This size cannot be directly measured from 
single epoch spectra of AGN, which are used for estimate 
of BH masses. Most popular estimated of the BLR size rely 
on the discovery that Rblr scales with a certain power of 
continuum luminosity of the AGN. Unfortunately this scal¬ 
ing dependence appears slightly different in works of various 
authors. 

For example, for Fairall 9 an d some othe r objec ts we 
used dat a for BH masses from iBrenneina 3 (120131 1 and 
[Reynold^ (l2013l l. that are based on the estimate of the 
virial coefficient in Eq.([T]) as f — 1.0. For Mrk 290 
iFeng. Shen fc Lil |(20l3) produced the estimates of its mass 
using the specific calibration of the coefficient of single epoch 
spectra in the local AGN sample. This circumstance may be 
the reason of the difference between our and others results. 

There is yet another way for determining the SMBH 
mass. Under the assumption that the motion of the gas in 
the BLR of AGNs is dominated by the gravitational influ¬ 
ence of the black hole, one can use the virial relation ([T]) to 
obtain Mbh- The radi us Rblr is determin ed usually with 
reverberation mapping dPeterson et al.ll2004 ) or i s estimated 
with the radius-luminosity relation I Bentzll2009ll . 

If one adopts a value for the virial coefficient of / = 1, 
the determined Mbh is commonly called as "the virial 
product” (VP). The results of ca lculations o f the virial 
product for AGN s are presented bv iGrier et al.l ll2013ll and 
IHo fc Klr3 (l2014l l. Eq.([3]) allows to obtain the following re¬ 
lation between the actual value of BH mass, virial product 
and sini: 


sin i = 


1 

2 


/ VP 
V Mbh / 


(5) 


Determining the value of sini from the polarimetric data 
one can deri ve the actua l value of Mbh- Us i ng da ta on VP 
published in iGrier et al.l ll2013t ): lHo fc Kiniil (l2014h , we esti¬ 


Table 3. The masses of SMBH in AGNs determined via measured 
polarization of broad H g emission and vir ial product data. Po- 
larizatio n data are f r om [ Smith et alJ ll2002li . Virial product data 
are from iHo fc KimI ||20 i1l 


Object 

Mbh / Mq 

Mbh / Mq 


(observations) 

(from Eq.JSj) 

Akn 120 

(3.23/:° |^) X 10® 

(2.6410 50) X 10® 

Fairall 9 

(1.34+0]®) X 10® 

(1.261]]^) X 10® 

Mrk 509 

(1.35 ±0.12) X 10® 

(1.26 ±0.0053) X 10® 

NGC 3516 

]^q8.06±0.27 

03 

NGC 3783 

]^q7.7±0.11 

7 .7+0.02 

20^'^ -0.04 

NGC 4051 

(1-64]0-67) X 10® 

(1.710 0) X 10® 

NGC 4593 

(8.25t®]®) X 10® 

(1.15l0;244) X 107 

NGC 5548 

(7.8410-53) X 10^ 

(1.4910101) X 10® 

NGC 6814 

(8.7llg) X 10® 

(6.410 01) X 10® 

NGC 7469 

( 3 . 47 II 45 ) X 10 '^ 

(4.251010) X 107 

PG 1211-1-143 

(2.1912-07) X 108 

(2.17l0;®4) X 10® 


mated the values of BH masses for AGNs presented in Table 
2. These results are presented in Table 3. They agree with 
the values in Table 2 within the uncertainties. 


4 CONCLUSIONS 

We presented a new method for estimating the SMBH 
mass. The method is based on the observations of polar¬ 
ization of broad fines emission, suggesting that the BLR 
is likely to be an optically thick fl a ttened configuration 
dKollatschnv. Zetzl fc Dietrichl l2006l : iKollatschnv fc Zetzll 
l2013l) . In this case the vir ial coefficient (Eq. ( Hi) c an be 
derived as / = 1/(4 sin^ i) dMcLure fc Dunlod 20 Pill . The 
value sini can be obtained from the polarization degree 
of radiation of optically thick plan e parallel atmosphere 
llChandrasekharlll950l : [Sobolevlll963ll . This result allows us 
to obtain an independent estimate for the central supermas- 
sive black hole mass from Eq. ©• The results of our calcula¬ 
tion of the masses of SMBH for specific AGN are presented 
in Table 2. Comparison of our results with the results of 
other authors, which used fixed values of the virial coeffi¬ 
cients, and frequently adopted the value of / = 1, shows 
significant differences for a number of AGNs. It means that 
there is a clear distinction of our black hole masses from 
the virial product value that correspond s to the value of the 
virial coefficient / = 1 iIHo fc Kimlf20l3 L 

It is important to emphasize that the polarimetric ob¬ 
servations allow to derive the ratio Rblr/R g for the situa¬ 
tion when the BLR resides in a flattened, disk like configu¬ 
ration (see Eq. ®)- 
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